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We demonstrate that an optically driven spin of a carrier in a quantum dot undergoes indirect 
dephasing via conditional optically induced charge evolution even in the absence of any direct 
interaction between the spin and its environment. A generic model for the indirect dephasing with 
a three-component system with spin, charge, and reservoir is proposed. This indirect decoherence 
channel is studied for the optical spin manipulation in a quantum dot with a microscopic description 
of the charge-phonon interaction taking into account its non-Markovian nature. 
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I. INTRODUCTION 

The idea of using quantum dots (QDs) for quantum 
computer implementations follows from the possibility 
of a clear selection of a two level system, on which a 
qubit can be realized To this end, both charge and 
spin states of confined carriers are employed, where the 
latter is preferable, since spin states are generally more 
resistant to decoherence processes. Moreover, it is pos- 
sible to exploit the charge evolution dependent on spin 
(via selection rules and Pauli exclusion principle) in or- 
der to manipulate the spin by optical means [2, H, 0, Q 
on picosecond time scales, that is, much faster than pre- 
viously proposed magnetic or electrical control. Many 
spin control schemes in such hybrid systems which use 
off-resonant interband excitations together with STIRAP 
processes Q , adiabatic 0] and fast [a, [gj] evolution within 
trapped states in A or four-level [l0[ systems have been 
proposed. These hybrid systems are considered now as 
the most promising candidates for QD-based quantum 
computers since during the millisecond spin decoherence 
time [11] it is possible to perform about 10 9 optical quan- 
tum gates. Optical rotation of a single spin performed 
via picosecond laser pulses with the optical Stark effect 
as the operative mechanism was recently experimentally 
demonstrated [12]. This pioneering experiment showed 
that fast optical spin control is feasible and the current 
task is to thoroughly study the decoherence mechanisms 
that limit the fidelity of the achieved quantum control. 
The fundamental question is whether the spin degrees of 
freedom are indeed affected by decoherence mechanisms 
to a smaller degree than the charge ones and what con- 
stitutes their main dephasing channel. 

In this paper, we show that the spin state of a con- 
fined carrier can undergo dephasing even in the absence 
of spin-reservoir coupling if the spin rotation is achieved 
by a conditional evolution induced on the orbital degrees 
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of freedom, as is the case in an optical control scheme. 
Although the dynamical details of this dephasing process 
depend on the specific implementation, the fundamental 
idea of the indirect dephasing can be understood with the 
help of a "generic model" of a three-component system: 
the carrier spin, its orbital state, and the reservoir. We 
show that this additional decoherence channel occurs on 
comparable or even shorter timescales than the spin pre- 
cession and trion decay during the optical manipulation. 
Thus, it may constitute the main source of imperfections 
of the optical spin rotations. This shows that phonon- 
induced dephasing should be included in the analysis of 
optical spin control schemes even though the commonly 
studied decoherence mechanism related to the material 
dependent spin-orbit coupling leads to very small errors 
for short gates [ll| and indeed can be neglected. 

The paper is organized as follows. In Sec. HI a generic 
model describing the indirect spin dephasing is intro- 
duced. Next, in Sec. IIII[ we present the model for the 
specific optical spin control protocol in a single QD. Sec- 
tion IIVI describes decoherence processes resulting from 
carrier-phonon coupling. Section [V] concludes the paper 
with final remarks. 



II. INDIRECT DEPHASING 

The idea of optical spin rotation is based on a spin- 
dependent evolution of the charge, which finally brings 
it to the original state, up to an additional phase accu- 
mulated during the evolution. Let the initial state be 
IV^o)) = (ck|0) s + f3\l) s ) |0) c , where the components 
refer to spin (s) and charge (c) states, respectively. The 
ideal evolution then has the form: 

|^id(*)> = a\0) B ® |0) c + p\l) s ® [f;(t)|0) c + £(t)|l)c] , 

where, at the final time ti, rj(ti) = e 1 ^ and = 0. 

Typically, the occupation of the excited charge state is 
kept small, <C \r]{t)\. This evolution realizes a rota- 

tion of the spin by an angle <fi around the axis defined by 
the states |0) s , |l) s , which may be selected at will using 
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selection rules and appropriate pulse phases and polar- 
izations. 

While the interaction between the spin and the envi- 
ronment is very weak, there is much stronger scattering of 
the reservoir quanta on the charge excitation, inducing, 
in a static situation, the usual phase damping channel 
on the charge subsystem. In the present case, when the 
charge state performs a conditional loop in its Hilbert 
space, the transient occupation of the excited charge 
state leads to the accumulated scattering amplitude (in 
the leading order in £ and e), 




where |e| 2 is proportional to the scattering rate and we 
assume that the reservoir quanta are non-resonant with 
the transitions between the charge states (otherwise, ad- 
ditional leakage out of the computational subspace ap- 
pears). The final state of the three-component system is 
therefore 

|^ac(tl)) = a|O) s 0|O) c 0|O) e 

+e^|l) s |0) c (Vl-M 2 |0)e + Hl)e), 

where the last component (e) represents the environment 
states. Thus, the charge state separates but the spin 
state becomes entangled with the environment. Tracing 
out the charge and environment degrees of freedom one 
arrives at the operator sum representation for the effect 
of the imperfect rotation on the spin state, 

l 

fi=0 

with M = |0) ss (0| + Vl-M 2 |l)ss(l|, Mi = M|l)ss(l|. 
In this way, the coupling between the orbital degrees 
of freedom and the reservoir has induced an indirect 
phase damping channel on the spin qubit (in the gate- 
dependent basis |Q) S , |l) s ), analogous to the indirect mea- 
surement scheme [13] with the spin, charge and environ- 
ment playing the roles of the quantum object, quantum 
probe and measurement device, respectively. 

In the following, we study in detail the indirect dephas- 
ing process for a specific optical spin control protocol Q , 
including the microscopic description of the interaction 
between charges and their phonon reservoir as well as 
the non-Markovian nature of the latter. We show that 
this dephasing process leads to considerable errors, much 
larger than those induced by the spin-orbit coupling or 
hyperfine interaction over the relatively short gate dura- 
tion on the picosecond timescale. 

III. MODEL SYSTEM 

The considered system consists of a single QD doped 
with one electron. A magnetic field is applied in the 
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FIG. 1: A system in a single quantum dot. 

x direction (Voigt configuration) and generates Zeeman 
splittings 2uo e between the two electron spin states \x) 
and \x) with fixed spin projection on the x axis equal to 
— 1/2 and +1/2, respectively. Analog for the trion spin 
states \Tx) and \T X ) with energy splitting 2cjh- These 
states are linear combinations of the electron (\z), \z)) 
and trion (|T), \T)) spin states along the growth and 
optical axis z. Depending on the light polarization, ro- 
tations about different axes are accomplished. 

As shown in Ref. |8|, a rotation about the z axis is per- 
formed with off-resonant circularly polarized light which, 
according to selection rules, couples the two spin states 
to only one trion state. Thus, we deal with an evolu- 
tion of a three- level A system (see Fig. [1]). The control 
Hamiltonian, including free carrier part and carrier-light 
interaction, reads 

H C = u; e (\z)(z\ + \z)(z\) + e T \T)(T\ 
+Q z (t) (e^ £ |z)(T|+H.c), 

where the laser pulse couples only the one spin state \z) 
and a trion state |T), whereas the orthogonal spin state 
\z) is indirectly coupled via the magnetic field. After a 
passage of a 2tt sech pulse, Q z (t) = Q z sech.(cr z t), the state 
acquires a phase, which, in consequence, leads to a spin 
rotation. The angle of rotation, <p z = 2 arctan(<j 2 /A 2 ), is 
defined via the laser bandwidth <j z and detuning of the 
laser from the transition energy A z = er — oo z . No pop- 
ulation transfer to a trion state is possible for <j z = Q z . 
The approximation made in this scheme requires that 
the spin is considered to be frozen during the pulse, i.e. 
cr z ^> 7, where 7 = 2(u e -f c^), which from the beginning 
imposes a limitation on driving conditions (short pulse 
durations especially for large Zeeman splittings). 

The free phonon Hamiltonian has the form H p h — 
^fc^fcAfcAfc> wnere is the phonon wave number and 
Ate (At) * s a ph° non creation (annihilation) operator with 
corresponding frequencies cok- The Hamiltonian describ- 
ing the interaction of the carriers with phonons reads 

ffc-ph = £ \n)(n'\ £/nn'(fe)(/Jfc +/£„), 
n,n' k 

where f nn >(k) are coupling elements and n = z, z, T, and 
T . The off-diagonal elements can be neglected due to 
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FIG. 2: Phonon spectral density R(co) at two temperatures 
and two spectral characteristics of the driving si(uj) and s 2(1*1) 
for 7r/2 rotation about the 2; axis with detuning and pulse 
bandwidth a z = A z = 2.6 meV. 



energetic reasons and low efficiency of direct phonon- 
assisted spin- flip processes. Moreover, f zz (k) = f zz (k) 
since the orbital wave functions are the same. Before 
the pulse is switched on, the lattice is already in a new 
dressed equilibrium state [14[ due to doping with one 
electron, and the phonon modes can be redefined in terms 
of new operators b k = flu + fzz(k)/(huok)- In the strong 
confinement regime, a trion state can be written in a 
product form of electron and hole states. The resulting 
carrier-phonon Hamiltonian is 

H c _ ph = \T)(T\J2F TT (k) (b k + bl k ) 

k 

with the following deformation potential coupling ele- 
ment between a trion and the phonon environment pj| 



F TT (k) = f TT (k) - f zz (k) 



I hk 
2pVci 



(D e -D h )T(k). 



Here, p = 5360 kg/m 3 is the crystal density, V is the nor- 
malization volume of the phonon modes, c\ = 5150 m/s 
is the longitudinal speed of sound, J~(k) is the form 
factor reflecting the geometrical properties of the wave 
functions [l6[, and D e (Dh) is the deformation poten- 
tial constant for electrons (holes), where D e — = 
8 eV. These parameters correspond to self-assembled 
InAs/GaAs quantum dots with the electron and hole con- 
finement in-plane equal to 4 nm and in growth direction 
1 nm. 



IV. PHONON-INDUCED DECOHERENCE 

To measure the quality of the operation on a qubit 
we use the error of the quantum gate, S = 1 — F 2 , 
defined as the loss of fidelity F. The error is a dif- 
ference between the ideal final state (without decoher- 
ence) and the actually achieved one including the cou- 
pling to environment. Here, we consider the interaction 
with phonon environment, however, the trion radiative 
coupling (carrier-photon interaction) can be described in 
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FIG. 3: Phonon-induced error contribution due to (a) pure 
dephasing and (b) phonon-assisted trion generation during 
the 7r/2 rotation about the z axis. 



the same manner. The effect of the interaction with the 
phonon reservoir is calculated via the second order Born 
expansion of the density matrix evolution equation (for 
details, see Ref. Il5h . The interaction with light is in- 
cluded exactly and coupling to phonons is treated within 
a non-Markovian perturbation theory. As a result, one 
can write the error of the quantum gate as an overlap 
between two spectral functions reflecting the properties 
of the two above interactions, 



= J oIujR(uj)S(uj). 



Here, 



n B 



h 2 



-E 



F TT (k)\ 2 [S(u; - cu k ) + 5{u + cj k )] 



is the phonon spectral density representing phonon emis- 
sion (lj > 0) and absorption (uo < 0, nonzero only at 
finite temperature) processes (see Fig. [2]). 

The spectral characteristics of the driving, S(u), has 
as many contributions as the dimension of the orthogonal 
complement of the initial state. In the case of z rotation, 
there are two contributions, S(co) = si(w) + 52 (00) re- 
flecting two phonon-induced decoherence channels. One 
represents pure dephasing mechanism and reads 



si(u) = - sin 2 $ 



I 



dt 



-iut 



--r (i 

c* 



■0' 



where c = (l + iA z /<j z )/2 and the time dependence is en- 
closed in £(t) = [tanh(cr^t) + l]/2. This function is always 
centered at uo = (Fig. [2]) and its width grows with grow- 
ing pulse bandwidth and detuning. It results from the 
fact, that the dynamical errors depend on the evolution 
speed, i.e., for a given pulse duration only some phonon 
modes can follow the evolution adiabatically whereas the 
others relax contributing to dephasing. The same applies 
to the second spectral function 
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cos — 
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dt e 



-iut 
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but the center of this function is shifted to a negative 
frequency around detuning value uo « — A z . This contri- 
bution represents real transition and constitutes a deco- 
herence channel referred to as the phonon-assisted trion 
generation. 
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FIG. 4: Total phonon-induced error for (a) positive and (b) 
negative detuning during the tt/2 rotation about the z axis. 
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FIG. 5: 4— level system in a single quantum dot. 



The resulting phonon-induced errors during a tt/2 rota- 
tion about the z axis, averaged over all initial spin states, 
are plotted in Fig. [3] as functions of detuning and band- 
width (in this case, A z = a z ) at four different tempera- 
tures T. 

The first contribution to the error resulting from pure 
dephasing effects [Fig. [31(a)] initially grows with growing 
detuning and pulse bandwidth. For small pulse band- 
widths, the evolution is really slow and the relevant func- 
tion Si(uj) is extremely narrow covering only the dimin- 
ishing part of the phonon density at oj « 0. Thus, the 
phonons are able to adiabatically follow the change of 
the charge distribution and as a result the decoherence 
is reduced. Unfortunately, the proposed schemes require 
usually bandwidths much larger than Zeeman splitting 
and one cannot use the discussed bandwidth sector with 
small errors. This error contribution reaches its maxi- 
mum value for A z = <j z w 1.5 meV for all temperature 
values, where the pure dephasing effects are most effi- 
cient [si(u) is broad and covers the whole spectrum of 
phonons] . 

The second error due to phonon-assisted transitions 
to the trion state is plotted in Fig. [3f b). In this case, 
the temperature dependence is stronger, since the spec- 
tral characteristics is centered at the negative frequency 
part of the phonon spectral density, which is strongly 
temperature dependent. Even for small bandwidths at a 
relatively low temperature T = 1 K, the error is larger 
than 10 -4 . At each temperature, the maximum error is 
reached for the detuning corresponding to the maximal 
value of the phonon density. The error diminishes for 
large detunings (> 50 meV) after the spectral character- 
istics reaches the phonon cut-off, where the one-phonon 
processes are not efficient. 

The total phonon-induced error during the tt/2 rota- 
tion about growth direction z is plotted in Fig. Ufa). 
To guarantee the coherent control and reach small er- 
rors, one needs either very small values of detunings and 
pulse bandwidth or very large ones of a few tens of meV. 
Taking into account the bandwidth limitation for typical 
Zeeman splitting of 0.1 meV, the available parameters 
lead to large gate errors even at zero temperature. The 
only way to obtain desired small errors is to use very 
large detunings and short pulse durations. However, un- 
der such conditions, many other decoherence channels 
like resonant and off-resonant transitions to higher states 



or interaction with optical phonons are likely to appear. 
Moreover, this can lead to experimental difficulties, since 
large detunings require very strong pulses. 

In order to perform a rotation about an arbitrary 
axis, rotations about two orthogonal axes are needed, 
e.g. z and x, and detunings above the energy gap may 
be needed. This leads to much larger phonon-induced 
errors, since emission processes become very important 
here. The total phonon-induced error for negative detun- 
ings is plotted in Fig. [H[b). Now, the spectral character- 
istics S2(oo) responsible for phonon-assisted trion genera- 
tion is centered at positive frequencies, where the phonon 
spectral density has much larger values especially at low 
temperature. In this case, the errors are up to two or- 
ders of magnitude larger in comparison with those for 
positive detunings. For experimentally reasonable values 
of detunings and pulse bandwidth, the error is always 
larger than 10 -2 and has the maximal value of « 10 _1 
for A z = <j z w 1 meV. 

The spin rotation about the x axis is realized via lin- 
early polarized tt x pulse. The relevant control Hamilto- 
nian is: 

Hex = ^(kXxl-lxXxO + e^lT.XT.I+e^lT.X^I 
+ [Q x (t)e-™°\\x)(T x \ + \x)(T x \) +H.c] , 

where = ct ± ujh- In this case, all four levels par- 
ticipate in the evolution (see Fig. [5]). In consequence, 
there are two paths for phonon-assisted trion generation 
with two different detunings. The resulting total phonon- 



10" 




0.01 0.1 
a x (meV) 



FIG. 6: Total phonon-induced error for the tt/2 rotation 
about the x axis. 
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induced gate error for the tt/2 rotation about the x axis 
is shown in Fig.[6j One can see, that already at T = 1 K, 
the error is always larger than 10 -4 and grows with grow- 
ing bandwidth. 

Adding the individual errors, it is possible to estimate 
the error of an arbitrary spin rotation. As we already 
discussed, even for a rotation about one of the axes, it 
is impossible to find driving conditions leading to errors 
smaller than 10 -4 , thus for the ^/-rotation the situation 
is even worse. Moreover, the calculated errors for single- 
qubit gates provide an estimation for the two-qubit spin 
gates employing, for instance, the electron hole exchange 
interaction in coupled QDs [171 ]. 

V. DISCUSSION AND CONCLUSIONS 

It has been shown that even in the absence of direct 
spin-reservoir coupling, the spin state of a confined car- 
rier is exposed to indirect dephasing through the entan- 
gling optically induced charge evolution. We have pro- 
posed a model for this indirect decoherence channel con- 
sisting of three components: spin, charge, and environ- 
ment. As an illustration, the optical spin manipulation 
in a single doped quantum dot has been considered. It 
was shown, that optical driving of such a system leads to 
a strong dynamical response of the lattice and to strong 
indirect dynamical phonon-induced decoherence channels 
for the spin degrees of freedom. 

Finally, we compare the considered optical spin control 
proposal with two previous schemes [y, [fl] . All of them 
use single or double quantum dots doped with one addi- 
tional electron and the excitation of intermediate trion 
state. One of them [6] makes use of stimulated Ra- 
man adiabatic passage (STIRAP) and is implemented 
in a double quantum dot. The main limitations of this 
proposal are slow adiabatic evolution requirement and 
necessity of electron transfer between two QDs and the 
delocalized hole state. The second one [7| is implemented 



in a single QD so that the two latter constrains are over- 
come. However, the evolution still has to be adiabatic. 
The proposal considered in this paper prevails over all 
the limitations discussed above, since the optical rota- 
tion is performed by means of fast laser pulses. However, 
this leads to larger phonon-induced errors 5 > 10 -3 even 
at low temperature, whereas the errors in the case of adi- 
abatic evolution [15|, LL8[ are at least one order of magni- 
tude smaller S < 10 -4 . On the other hand, the fast evolu- 
tion leads to smaller errors resulting from carrier-photon 
interaction. The trion state is excited only for a short 
moment, thus the probability of its radiative decay is low 
[15[ . All in all, the fast optical spin rotation analyzed here 
possesses many advantages in comparison with the other 
two proposals, however the dynamical phonon-induced 
indirect spin dephasing is in this case much stronger. 

The phonon-induced decoherence processes may in 
many cases constitute the dominant source of errors, 
since they are much more efficient than those due to spin- 
orbit mechanism assisted by phonons and up to two or- 
ders of magnitude larger than the errors resulting from 
trion radiative decay [15]. Moreover, these dynamical 
phonon-induced processes are most efficient exactly on 
the timescales for the proposed and demonstrated opti- 
cal spin rotations. Therefore, in order to overcome the 
phonon-induced indirect spin dephasing one should avoid 
such detunings and timescales. Another idea is to reduce 
the dephasing by means of collective encoding of quan- 
tum information in QD arrays pjl l2(ill. Puls e op timiza- 
tion may also lead to error reduction [2l], [22|, [23| • 
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